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Abstract 

The problem of Bloch electrons in two dimensions subject to mag- 
netic and intense electric fields is investigated. Magnetic translations, 
electric evolution and energy translation operators are used to spec- 
ify the solutions of the Schrodinger equation. For rational values of 
the magnetic flux quanta per unit cell and commensurate orientations 
of the electric field relative to the original lattice, an extended su- 
perlattice can be defined and a complete set of mutually commuting 
space-time symmetry operators is obtained. Dynamics of the system 
is governed by a finite difference equation that exactly includes the 
effects of: an arbitrary periodic potential, an electric field orientated 
in a commensurable direction of the lattice, and coupling between 
Landau levels. A weak periodic potential broadens each Landau level 
in a series of minibands, separated by the corresponding minigaps. 
The addition of the electric field induces a series of avoided and exact 
crossing of the quasienergies, for sufficiently strong electric field the 
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spectrum evolves into equally spaced discreet levels, in this "magnetic 
Stark ladder" the energy separation is an integer multiple of hE/aB, 
with a the lattice parameter. 

1 Introduction 

The problem of electrons moving under the simultaneous influence of a pe- 
riodic potential and a magnetic field has been discussed by many authors 
[U[2Hnil31EHnHZllH] : the spectrum displays an amazing complexity includ- 
ing various kinds of scaling and a Cantor set structure [§]. Some of these 
results had been used in order to disclose the topological structure of the 
Hall conductance within the linear response theory fUJ El IE] • The addi- 
tion of an electric field leads to interesting new phenomena that makes its 
analysis worthwhile. For example, it has been suggested that as the strength 
of the electric field increases the longitudinal quasi-momentum is quantized 
leading to the appearance of a "magnetic Stark ladder" , in which the bands 
are replaced by a series of quasi discreet levels [13] . The quantization of the 
spectrum also affects the the localization properties of electron wave packets 

[Illlig. 

We consider the problem of an electron moving in a two-dimensional lat- 
tice in the presence of applied magnetic and electric fields. We refer to this 
as the electric-magnetic Bloch problem (EMB). The corresponding magnetic 
Bloch system (MB) has a long and rich history. An early important contri- 
bution was made by Peierls [I] who suggested the substitution of the Bloch 
index k by the operator (p — eA) in the B = dispersion relation e(fe), which 
is then treated as an effective single-band Hamiltonian. The symmetries of 
the MB problem were analyzed by Zak [SI El E], who worked out the repre- 
sentation theory of the magnetic-translations group. The renowned Harper 
equation was derived assuming a tight-binding approximation in which the 
magnetic field acts as a perturbation that splits the Bloch bands [2]. Rauh 
derived a dual Harper equation in the opposite limit of intense magnetic 
field [S] , here the periodic potential acts as a perturbation that broadens the 
Landau levels, in this case the Harper equation takes the form 

c m _i + 2 cos (2-Kam + n)c m + c m+1 = ec m . (1) 

where a = 1/0 is the inverse of magnetic flux in a cell in units of h/e. 
The studies of the butterfly spectrum by Hofstadter and others have since 
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created an unceasing interest in the problem because of the beautiful self- 
similar structure of the butterfly spectrum [HJ . Remarkably, an experi- 
mental realization of the Hofstadter butterfly was achieved considering the 
transmission of microwaves ^7] and acoustic waves ^B] through an array of 
scatterers. However, the study of these phenomena in electronic transport 
has only recently become experimentally accessible with the developments 
in the fabrication of antidot arrays in lateral superlattices by ion beam and 
atomic force lithography ^JJ EOj- For example, a signal of the Hofstadter 
Butterfly spectrum through the measurement of the magnetoresistance and 
Hall conductance in artificial arrays of anti quantum dots has just been re- 
ported [2*T] . It is also worthwhile to point out that Aubry and Andre f22] 
utilized Eq. [T] to model a one dimensional system that develops a metal- 
insulator transition; this transition, arises from the competing effects of two 
potentials terms, which are periodic but incommensurate with one another. 

The symmetries of the EMB problem were analyzed by Ashby and Miller, 
who constructed the group of the electric-magnetic translation operators, 
and worked out their irreducible representations 23 . The properties of the 
electric-magnetic operators were utilized in order to derive a finite difference 
equation that governs the dynamics of the EMB problem when the coupling 
between Landau levels can be neglected [15] . In this paper we shall derive 
the equation that applies under most general conditions. Magnetic transla- 
tions, electric evolution and energy translation operators are used to specify 
the solutions of the Schrodinger equation, commensurability conditions must 
be implemented in order to obtain a set of mutually commuting space-time 
symmetry operators. In addition to the broadening of the Landau levels pro- 
duced by the periodic potential, the electric field induces a series of avoided 
and exact level crossing that yields a generalized Hofstadter spectrum, with 
an interesting pattern that arises out of the competition between by the lat- 
tice potential and the external E and B fields. For sufficiently strong electric 
field the spectrum evolves into equally spaced discreet levels, in this "mag- 
netic Stark ladder" the energy separation is an integer multiple of hE/aB, 
with a the lattice parameter. 

The paper is organized as follows. In the next section we present the 
model that describes the EMB problem and construct its symmetry opera- 
tors. In Section |3] we describe the commensurability conditions required to 
have simultaneously commuting operators, they are exploited in order to con- 
struct an appropriated wave function basis. We derive an effective equation 
(|3U|) . in which the "evolution" is determined by a differential equation with 
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respect to the longitudinal pseudomomentum. The derivation of the finite 
difference equation that governs the dynamics of the system is presented in 
Section HJ Results for the energy spectrum are presented, and are also dis- 
cussed from the perspective of the adiabatic approximation. The last section 
contains a summary of our main results. 



2 The electric-magnetic Bloch problem 
2.1 The model 



Let us consider the motion of an electron in a two-dimensional periodic 
potential V, subject to a uniform magnetic field B perpendicular to the 
plane and to a constant electric field E, lying on the plane according to 
E = E(cos6, sin 8) with 9 the angle between E and the lattice x— axis. The 
dynamics of the electron is governed by a time-dependent Schrodinger equa- 
tion that for convenience is written as 



here m* is the effective electron mass, IT^ = p^+eA^, with p^ = (ihd/ dt, —iKV) 
Throughout the paper energy and lengths are measured in units of fiw c = 



and £ = \l jk, respectively, where uj c is the cyclotron frequency and £q the 



magnetic length. So, unless specified, we set h = e = m* = 1; although B 
can also be omitted from the expressions, we find convenient to explicitly 
display it. Covariant notation will be used to simplify the expressions, e.g. 
%n = (t,x) = (t,x,y). 

Equation (J2J) can be considered as an eigenvalue equation for the operator 
S with eigenvalue 0. The gauge potential is written in an arbitrary gauge, 
it includes two gauge parameters a, and /5, the final physical results should 
be of course, independent of the gauge. Hence the components of the gauge 



s\1>)= ^(^ + n;)+v-n |^) = o, 



(2) 
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potentials are written as 



A = 
A™ — 



Ay — 



P-^) x-E, 



P+^E x t+(a-^)By, 



(3) 



The symmetrical gauge is recovered for a = 0, whereas the Landau gauge 
corresponds to the selection a = 1/2. If /? = 1/2 all the electric field contri- 
bution appears in the vector potential, instead if (3 = —1/2 it lies in the scalar 
potential and the Schrodinger equation becomes time independent, however 
even in this case a time dependence will slip into the problem through the 
symmetry operators. A general two-dimensional periodic potential can be 
represented in terms of its Fourier decomposition 



v i x , V) = v ™ ex P ( 1 1" % J 

r,s ^ ' 



(4) 



For specific numerical results we shall use the potential 

V(x,y) = U 



27rx\ ( 2ny 

cos | I + A cos 



(5) 

V u / 

A is a parameter that can be varied in order to have an anisotropic lattice; 
A = 1 corresponds to the isotropic limit. 



2.2 Electric evolution and magnetic translations 

Let (t, x) — > (t + r, x + R) be a uniform translation in space and time, 
where r is an arbitrary time and R is a lattice vector. The classical equa- 
tions of motion remain invariant under these transformations; whereas the 
Schrodinger equation does not, the reason being the space and time depen- 
dence of the gauge potentials. Nevertheless, quantum dynamics of the system 
remain invariant under the combined action of space-time translations and 
gauge transformations. The electric and magnetic translation operators are 
defined as 

T (r) = exp (-itO q ) , T,j(a) = exp (ia Oj) , (6) 
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with j = x,y and the electric-magnetic symmetry generators are written as 
covariant derivatives M = + A M , with the components of the dual gauge 
potentials A M given by 

A = A + x ■ E , 

A x = A x + By + E x t , A y = A y - Bx + E y t . (7) 

It is straightforward to prove that the operators in (JSJ are indeed symmetries 
of the Schrodinger equation; they commute with the operator 5* in Eq. @. 
Similar expressions for the electric-magnetic operators were given by Ashby 
and Miller [23], however their definition included simultaneous space and 
time translations; we deemed it more convenient to separate the effect of the 
time evolution generated by the To to that of the space translations generated 
by Tj. The following commutators can be worked out 

[Uo,^] =-iEj, \n u U 2 ] = -iB, 

[0 Q ,O j ] = iE j , [O u 2 ]=iB, 

[n M ,a] = o. (8) 

Notice that the electric-magnetic generators have been defined in such a 
form that they commute with all the velocity operators IT. The Schrodinger 
equation and the symmetry operators are expressed in terms of covariant 
derivatives 11^ and O^, respectively. A dual situation in which the roles of 11^ 
and Op, are interchanged can be considered. The dual problem corresponds 
to a simultaneous reverse in the directions of B and E. We notice that 
the commutators in the second line of Eq. (jHj) are part of the Lie algebra 
of the EM-Galilean two dimensional group j2U This group is obtained 
when the usual rotation and boost operators of the planar-Galilean group 
are replaced by their electric-magnetic generalization in which the operators 
are enlarged by the effect of a gauge transformation. 

3 Electric-magnetic Bloch Functions 

3.1 Commensurability conditions 

In order to construct a complete base that expands the wave function we 
require the symmetry operators to commute with each other. However we 
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have 



T,T v = e- c ^°^T v T^ (9) 

where c° = — r, and c 1 = c 2 = a. A set of simultaneously commuting 
symmetry operator can be found if appropriated commensurate conditions 
are imposed, we follow a three step method to find them: 

1. First we consider a frame rotated at angle 9, with axis along the lon- 
gitudinal and transverse direction relative to the electric field. An 
orthonormal basis for this frame is given by = (cos9,sin9), = 
(—sin9, cos9) and e 3 = e L x e T . The electric field is parallel to e L and 
the magnetic field points along e 3 . We assume a particular orientation 
of the electric field, for which the following condition holds 

tan^ = |^ = ^, (10) 
E x mi' K ' 

where m\ and m 2 are relatively prime integers. This condition insures 
that spatial periodicity is also found both along the transverse and the 
longitudinal directions. Hence, we define a rotated lattice spanned by 
the longitudinal = be^ and transverse 6<r — be? vectors, where b = 
ay/ ml + mf. The spatial components of the symmetry generator O are 
projected along the longitudinal and transverse directions: Ol = ^l-O 
and O t = e T ■ O. It is readily verified that [O , O T }=0. 

2. For the rotated lattice, we regard the number of flux quanta per unit 
cell to be a rational number p/q, that is 

# S I = ^ = E. (ii) 

We can then define a extended superlattice. A rectangle made of q ad- 
jacent lattice cells of side b contains an integer number of flux quanta. 
The basis vectors of the superlattice are chosen as qb^ and bx- Under 
these conditions the longitudinal and transverse magnetic translations 
Tjj (qb) = exp [iqbOj)j and T T (b) = exp (ibOr) define commuting sym- 
metries under displacements qb L and b T . Henceforth we shall either use 
the subindex (L,T) or (1,2) to label the longitudinal and transverse 
directions. 
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3. We observe that T and T L (qb) commute with T T . Yet they fail to 
commute with each other: 



T (r) T L (qb) = e**" E T L (qb) T (r) . (12) 

However the operators To and Ti(qb) will commute with one another 
by restricting time, in the evolution operator, to discreet values with 
period 

2tt 1 / b \ 

T = nr Q , neZ, r = — — = - — , (13) 

qbE p \v D ) 

where the drift velocity is vd = E/B and we utilized Eq. (JTTjl to write 
the second equality, b/vn is the period of time it takes an electron 
with drift velocity Vp to travel between lattice points. Eq. (JT3J) can 
be interpreted as the condition that the ratio of two energy scales is 
an integer: as we shall see (abE) is the magnetic-Stark ladder spacing, 
whereas (27rt>£>/&) is the quasienergy Brillouin width; hence Eq. (JT3J) 
represents the ratio of these two quantities. 

We henceforth consider that the three conditions (fTUj). (fTTj) . and (fT3|) hold 
simultaneously. In this case the three EM operators: the electric evolution 
Tq = T (t ), and the magnetic translations T L = T L (qb), and T T = T T (b) 
form a set of mutually commuting symmetry operators. In addition to the 
symmetry operators it is convenient to define the energy translation operator 

m 

Te = exp (-iy^J > ( 14 ) 

that produces a finite translation in energy by 27r/r = qbE. T E commutes 
with the three symmetry operators but not with S. Its eigenfunctions 

T E ?p = e iqbm ip, (15) 

define a quasitime d modulo r . 

3.2 Wave function and generalized Bloch conditions 

Having defined T , 72 and 7p that commute with each other and also with 
S, it is possible to seek for solutions of the Schrodinger equation labeled by 
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the quasienergy {£) and the longitudinal (ki) and transverse (k 2 ) quasimo- 
mentum according to 



T \£,k 1 ,k 2 ) = e' lT ° £ \£,k 1 ,k 2 ), 

T L \£,k 1 ,k 2 ) = e^\£,k 1 ,k 2 ), (16) 
T T \£, k h k 2 ) = e ik2b \£, h, k 2 ) . 

The magnetic Brillouin zone (MBZ) is defined by k\ G [0, 27r /qb] and k 2 e 
[0, 2-jt/b]. Similarly the quasienergy £ is defined modulo 27r/r = qbE. If a 
restricted energy scheme is selected for the energy, the first energy Brillouin 
region is defined by the condition £ e [0, 2tt/tq]. We shall find convenient to 
perform a canonical transformation to new variables according to 

Qo = t, P Q = O - l -E\ 

Qi = n 2 + p, Pi = n x , (17) 

Q 2 = 1 - Et, P 2 = 2 + E, 

that satisfy the commutation rules P u ] = —ig^ vi g^ u = Diag(—1, 1,1). 
Applied to Eq. (J2J) the transformation yields for the Schrodinger equation 



Hlib) , H 



±(p? + Ql)+V(x,y)-EP 2 



In this equation (x, y) must be expressed in terms of the new variables: 

a y = ^ {Ql -P 2) + f iQ2 - Pl) . (19) 

a b b 

On the other hand, the EM-symmetry operators and the energy translation 
operator take the form: 

T = e- lToPo , T L = e ^(Q2+BQo) ; (20) 



T T = e ibP \ T P =e~^ Qo 



Notice that these operators do not depend on the variables Pi, Qi- Thus, it is 
natural to split the phase space (Q^, P^) in the (Qi, Pi) and the (Qo, Po] Q2, P2) 
variables. For the first set a harmonic oscillator base is used. Whereas for 
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the subspace generated by the variables (Q , P ; Q 2 , P 2 ), we observe that the 
operators (7^, T[, Tj? , T E ) with all possible integer values of k, I) form a 
complete set of operators. The demonstration follows similar steps as those 
presented by Zak in reference j^Zj- Hence a complete set of functions, for the 
subspace (Qo, Pq] Q 2 , P2), is provided by the eigenfunctions of the operators 
(7q, T[, Trp, T E ). As starting point, we select a base of eigenvalues of the 
following operators 

A ] A \n, £, k 2 ) =n\n,£ 1 k 2 ) , 
P £, k 2 ) = £ \i-i,£ , k 2 ) , (21) 
P 2 £, k 2 ) = k 2 £, k 2 ) , 

where A and A^ are the lowering and raising operators of the /i-Landau level: 



A=-j=(P 1 -iQ 1 ), A^ = -^=(P 1 + iQ 1 ). (22) 

It is straightforward to verify that these states fulfill the required eigenfunc- 
tion condition (|16j) for the % and Tt operators. On the other hand Tl 
induces a shift in the £ and k 2 labels, while Te induces a shift qbE in the £ 
label. The previous considerations suggest that a state £, ki, k 2 ) that is 
a simultaneous eigenfunction of the four operators T E , %, T L and T T with 
the required eigenvalues (Eqs. ITB1 and ITBj) . can be constructed as a linear 
superposition of states of the form [T L ] 1 [T^T^ 1 ]" 1 \(J,,£, k 2 ). We write down 
the state, and verify their correctness: 

|A £, h, k 2 ) = [r L e- iqbkl ] 1 <n [e^e^ m - k ^] m \n, £, k 2 ) . (23) 

l fj,,m 

It is easy to check that this function satisfies the three eigenvalue equations 
for the symmetry operators in (|16p. Additionally it can be verified that the 
eigenvalue condition ()15j) for the energy translation operator is also fulfilled 
by imposing the periodicity condition c^ +p = c^. Every state in (|2*3^1 yields 
a set of different eigenvalues, a condition that follows from the fact that the 
selected set of operators is complete, hence the base of eigenvalues in (f2~H|) is 
complete and orthonormal. 

Up to this point we have constructed a base for the set of four (T , 72, 
Tt, T~e) operators, however the operator Te is not a symmetry of the problem, 
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consequently the solution of the Schrodinger is constructed as a superposition 
of all the states characterized by the state becomes 

\£, h, k 2 ) = I ' ddc (0) |a e, h, k 2 ) = [r L e- iqbki ] 1 b^ b ^- k ^ m \n, £, h) 

I fi,m 

(24) 

where 

b^= [ d$C (tf) c^e iqbmm . (25) 



Equation (j2T|) represents the correct state that satisfies the eigenvalue equa- 
tions (|16j). It is convenient to recast it in a compact form as 

\£,k) = W(k 1 )\£,k 2 ), (26) 

where \£,k) = \£,ki,k 2 ), and the operator W (ki) is defined as 

W (ifei) = ^ pie - * 5 ** 1 ]' = J2e iqbl{EQo+Q2 - kl) . (27) 



I' 

z i 
whereas the ket \£, k 2 ) is given by 



\E, k 2 ) = e ir7bm{Q ^b^ | M) £, k 2 ) . (28) 



fi,m 



The previous expression stresses the fact that the quantity e~ lcrbmkl bf^ does 
not depends on k±, this will be demonstrated below Eq. (jJSJ. These results 
will prove to be very useful to analyze the EM problem, in particular they 
allow us to obtain an effective Schrodinger equation where the "dynamics" 
is governed by the derivative with respect to the longitudinal pseudomomen- 
tum. Utilizing Eqs. (|26|) y ([27]) it is straightforward to prove the following 
relation 

P \£, h, k 2 ) = [P , W] \£, k 2 )+WP \£, k 2 ) = {£-iE- V fc ) \£, h, k 2 ) . (29) 
Thus the Schrodinger (fTH|) can be recast as 

(£-iE- V fc ) \£, h, k 2 ) = H \£, h, k 2 ) . (30) 
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As it will be lately discussed, this form of the Schrodinger equation becomes 
very useful to prove various properties of the system, in particular the adia- 
batic solution of the problem. 

The wave function takes a simply form if we adopt the (Pq, Pi, P%) repre- 
sentation. In this case ^fc (p) = (Pq, Pi, P 2 \S, k\, k 2 ) yields 



P 2 -k 2 + T (m + pl) 
(31) 

where M (Pi) is the harmonic oscillator function in the Pi representation 
M (Pi) = (Pi H = J_ = e- p2 ^H, (Pi) , (32) 

and (Pi) is the Hermite polynomial. 

On the other hand, it is sometimes useful to restore the space-time rep- 
resentation of the wave function. This problem provides a good example of 
the use of canonical transformations in quantum mechanics. Identifying the 
matrices that relate the (x^p^) variables to (Q^,P^j, it is possible to apply 
the method of reference [2H] hi order to obtain the desired transformation. 
The final result for the wave function in the space-time representation yields 

^ fc (t, x) = e ik - x - m u k (t, x) , (33) 
where the modulation function u(t, x) is given by 

Uk U^ x \ = 1 c -ixi\(a-l/2)x 2 +ki] c i(P+l/2)Etxi 

v 2tt 



Y^i^ m e- %qbElt e^- x ^ m+pl ^, (x! -k 2 + 2vr^±^ 

(34) 



x 

film 



here (p^ is the same harmonic oscillator function Eq. (j32j) . but now evaluated 
in the space representation. It is worthwhile to notice, that whereas there is 
not explicit gauge dependence for the wave function in the P representation 
(Eq. (jMJ)), the dependence on the gauge parameters a and (3 becomes explicit 
in the space-time representation (Eq. ()34j)). 
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Solution (|33I34|) includes a superposition of Landau-type solutions orig- 
inated beneath the spatial and time periodicity. The spatial periodicity is 
simply related to the external potential V. Whereas time periodicity arises 
from the conditions imposed to the symmetry operators in order to produce 
commuting symmetries; as discussed in Eq. (JT3j), the period is given by the 
time that takes an electron to drift between contiguous lattice points. Notice 
that Eqs. (|33I34|) follows from the Bloch and Floquet theorem. However in 
the electric-magnetic case the modulation functions u (t, x) are not strictly 
periodic, instead they satisfy the generalized Bloch conditions 

u(t + t , xi,x 2 ) = e lT ° A ° u(t, xi, x 2 ), 

u(t, xi + qb, x 2 ) = e' iqbAl u(t, x 1: x 2 ), (35) 

u(t, xi, x 2 + b) = e ~ lbA ' 2 u(t, xi,x 2 ), 

the phases are determined by the dual-gauge potential that appears in the 
symmetry operators Eq. (J7j). It is straightforward to verify that Uk in 
Eq. (|34j) satisfies these conditions. The generalized Bloch conditions for 
the magnetic-Bloch problem have been previously discussed, our results in 
(|33j) reduces to those in reference |7], when E = and the gauge a = 1/2 is 
selected. Our results in (|3*3j) have been extended in order to include the elec- 
tric field effects, furthermore they apply for an arbitrary gauge. The correct 
normalization conditions for the modulation function are obtained as follows 



(27T) 2 



d x u* (£, x) u (t, x) — 1 , (36) 

MUC 



' J2 

qb 2 

where MUC represents the magnetic cell defined by: X\ G [0, qb] and x 2 G 
[0,6]. 

The function u satisfies conditions similar to those in Eq. (J35|) but in the 
MBZ. The MBZ is actually a torus T 2 , so the edges (ki = 0, k 2 ) and (fci = 
27r/qb, k 2 ) must be identified as the same set of points, the wave function can 
differ at most by a total phase factor (similarly for the edges (ki, k 2 = 0) and 
{k 1: k 2 = 2<K/b)): 

u(k 1 + 2ir/qb,k 2 ) = e %h u{k x , k 2 ), u(k t , k 2 + 2-n/b) = t %h u{k x , k 2 ), 

(37) 

the functions f\ (k) and f 2 (k) can be related with the quantized values of 
the Hall conductance ^2], however instead of Eq. (|33j) . these functions are 
not analytically known and must be computed numerically. 
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4 Harper generalized equation 



4.1 Finite difference equation 



The coefficient bf^ in Eq. (j24j) satisfies a recurrence relation that is obtained 
when this base is used to calculate the matrix elements of the Schrodinger 
equation (|18|). First, let us consider the contribution arising from the periodic 
potential in Eq. The x and y coordinates are written in terms of the new 
variables (Qi, P\,Q2, P2) by means of Eq. (|19j). producing a term of the form 
exp [(i27irrrii(Q 1 — P 2 )/b], and another contribution in which Qi and P 2 are 
replaced by Q 2 and P\. Once that Qi and Pi are replaced by the raising and 
lowering operators in Eq. (|22|). one is lead to evaluate the matrix elements 
of the operator D = exp (zA^ — z*A) that generates coherent Landau states. 
A calculation yields 

(2) = (z/ |exp (zA f - zM) I n) = 



%I4T" (\z\ 2 ) , //>^, 

S^r" (n 2 ) > a* < ^ 

(38) 

where are the generalized Laguerre polynomial. The evaluation of the 
terms that include the P 2 operator is direct, because the base is a eigenvalue 
of this operator, on the other hand the operator exp [(i2irrmiQ2) /b] acts 
as a translation operator that produces a shift on the index m of the bf^ 
coefficient. Taking into account these results, it is possible to demonstrate 
after a lengthly calculation that the Schrodinger equation (fTHj) becomes 

t=N 

A f m (ki, k 2 ) b m+t = {£ + Ek 2 + obEm) ~b m . (39) 

t=—N 

Besides its dependence on the index m, b m is a vector with L components 
and (ki, k 2 ) is a L x L matrix in the Landau space, according to 

bm — {b m } = ip m , b m) 6 m ) , 
"*)"" = e- Mt B % (r, s) 5 t , m2 . smi . (40) 



in 



In the previous expressions L is the highest Landau level included in the 
calculations, iV = max{(r, s) (mi + m 2 )} is related to the largest harmonic 
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in the Fourier decomposition of the periodic potential (J1J, and 
given by 



m 



r, s) 



{v 00 + H + |) 5^, r, s = 0, 

V rs D l>IX (H rs ) Q iK rs e iM rs [ab(rm2-srni+-m)+k2] r g _^ Q 

the following definitions where introduced in the previous equation 
V2ti 



(41) 



K r 



b 

2tt2 

If 
2tt 



(m 2 — mi) (ir + s) , 
(rm.2 — sroi) (rmi + sm2) 



(42) 



(roi! + sm 2 ) 



We recall that the indexes (r, /) refer to the Fourier expansion of the pe- 
riodic potential Eq. (@J), whereas (m 1 ,m 2 ) correspond to the integers that 
determined the electric field orientation in Eq. (fTU|) . 

Equation ()39|) describes the system in terms of a recurrence relation on 
a two dimensional (fi, m) basis. The coupling in m ranges from m — N to 
m + N. We follow the method of reference to recast the recurrence 
relation on m with N nearest-neighbor coupling into a tridiagonal vector 
recurrence relation. This is enforced by defining a N component vector 



C m — ( &jVm-lj fr/Vrnj ^Nm+li 

and matrices Q^j and with elements 



? iV(m+l)-l 



) 



(43) 



/ A"" 



Q m {ki,k 2 ) 



Nm 





v -iV+l 
-N 



A-l 



x iV(m+l)-l / 



<5m (fcl, fe) 



A 

Nm-\-l 



A 1 

AO 



aJV-1 

A^-2 
JVm+1 



. Al-JV A2-JV 
\ **iV(m+l)-l **N(m+l)-l 



A 



(44) 



iV(m+l)-l / 
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Using the fact that Q+ = (Q~V, Eq. flU can be reorganized tridiagonal 
recurrence relation 

Q~ c m _i + Q m c m + Q^m+i = [(£ + Ek 2 + abEm) I N + abED N ] c m . (45) 

Here, c m is iV— dimensional vector according to the definition (}4"3*|) . addi- 
tionally each one of its components is a L dimensional vector due to its 
dependence on the Landau index Eq. (|4L)p. I N is the unit matrix and 
Dn = Diag (0, 1, 2, . . . , N — 1). Likewise, each of the components of the 
matrices Q m , In, and Dn are L x L matrices relative to the Landau con- 
tributions. In relation with Eq. ()28j) it was mentioned that the quantity 
e -iabk x does not explicitly depend on the longitudinal pseudomomentum 
component ki, this can be easily verified if the substitution c' m = e~ tabkl c m 
is incorporated into Eq. (jMJ), then using Eqs. fljOj) and (jSJ), it is straight- 
forward to verify that is indeed independent of k±. 

Eq. is one of our main results, and deserves to be emphasized. It is 
a generalization of the Harper equation that exactly includes the following 
effects: (1) an arbitrary periodic potential, (2) an electric field orientated in a 
commensurable direction of the lattice, and (3) the coupling between different 
Landau levels. So far, no approximations are involved, consequently it holds 
under most general conditions. In practice, only a finite number L of Landau 
levels can be included on the calculations, but a very good convergence can 
be obtained with a reasonable small selection for L. Previous known results 
are recovered if some approximations are enforced: (i) if the electric field is 
switched off, its orientation is meaningless, so the integers (m 1; m 2 ) can be set 
to m% = 1 and m 2 = 0, in this case Eq. (J45j) reduces to the model previously 
discussed by Petschel and Geisel jSU]. If E = and additionally the 
coupling between different Landau bands is neglected and the potential is 
taken as the sum of cosines given in Eq. (j3J), then the system reduces to a 
set of Harper equations ((TJ), one for each Landau band. 

Let us return to the case that includes the electric field. Equation (|43|) 
describes the system dynamics under most generals conditions; in spite of its 
complicated structure the equation can be solved by the matrix continued 
fraction methods. However, a considerable simplification is obtained if the 
electric field is aligned along the a;— lattice axis, i.e. (mi = 1, 777,2 = 0). 
Henceforth, we consider that the electric field is orientated along the x— axis 
(mi — 1; m 2 — 0), so b = a and additionally that the periodic potential takes 
the form given in (j3J). The dimension iV of the vectors c m and matrices in 
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Eqs. (143)) and ()44|) reduces to N = max{(r, s)(mi + m 2 } = 1, hence the 
dimensions of c m = c m and Q reduce to include only the Landau indexes. 
Using equations (j4*Uj) . (j4"T]) and (j4*4"j) the matrices Q become 

(Q+r = (Q + r = 2{ il K x)a 2 etabklD " u (v 7 ^) , 

ttK 

{QmT = on , rs 2 e i(27r< ™ (0^) + cc, (46) 

Z ^1 -\- \j CL 

here the parameter is a measure of the strength of the coupling between 
Landau bands 

ma 2 U (1 + A) 

x = — — • (47) 

Taking into account these simplifications and that cancels out, we find 
that Eq. (|4"5j) reduces to 

<5"c m _i + QmCm + Q + c m+1 = (S + Ek 2 + abEm) c m . (48) 



4.2 Numerical results 

The generalized Harper equation is given by a tridiagonal infinite recurrence 
relation (|48|). If .E is switched off, the equation becomes periodic with period 
p, in this case the equation can be recast as a finite pL x pL matrix that can 
be solved by a direct diagonalization. However, the introduction of E breaks 
the periodicity, and most of the methods used to solve the Harper equation 
break down. Eq. (|48|) was solved using an expansion of the associated Green's 
function into matrix continued fractions (MCFs) [29J. The energy spectrum 
is determined detecting the change of sign in the Green's functions that 
appear in the vicinity of a pole. The density of states can be obtained from 
the expression 

N (S) = -Tr G {£) = T-Im [Tr G ± {£)} , (49) 

7T 7T 

where the discontinuity in Green's function is given by G (£) = G + (£) — 
G~ (£), and the retarded and advanced Green's functions are defined on the 
side limits G ± {£) = lim e ^ G {£ ± ie). The numerical solution is obtained by 
truncating the iteration of the (MCFs) after the M— th term. The solution 
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converges if M is large enough, in the calculations we observe that in order 
to obtain a convergence with a precision of one part in 10 7 , the cutoff can be 
selected as M ^ 100. 

In our calculations we have used the effective mass m* = 0.067m e typical 
for electrons in GaAs and a superlattice constant a = 100 nm [21]. The 
rescaled energy spectrum 8/ [L^e -7 ™"/ 2 ] is shown in Fig. Q for the lowest 
Landau level, a weak modulation is considered (JJq = 0.5meV) so the p — v 
Landau mixing is negligible. The electric field intensity is E = 0.5V/ cm cor- 
responding to a ratio of the electric to the periodic potential of p = eaE/U = 
0.02. In the strong magnetic region, o G [0,1], the Hofstadter butterfly is 
clearly depicted. A distorted replica of the butterfly spectrum can be still 
observed in the region u g [1,2]. As the magnetic field loses intensity, the 
effect of the electric field becomes dominant, the butterfly is replaced by 
discreet levels separated in regular intervals, i.e. a Stark ladder. Although, 
p = eaE/Uo is small, there is a regime in which the electric field dominates 
over the lattice potential contribution, we can trace down the origin of this 
effect to Eqs. flHBj) and showing that the periodic contribution is modu- 
lated by a factor e -7r(T ' 2 . Consequently, in addition to a small Landau mixing, 
the condition to preserve the butterfly spectrum can be stated as 

eaE < U e- na/2 . (50) 

This condition restricts the intensity of both the electric and magnetic fields, 
these are estimated as (inserting units): E <C U /(ea) ~ 100 (U /meV) V /cm 
and B > nh/(ea 2 ) ~ 0.6T. 

Fig. |21 shows the energy spectrum for the four lowest Landau bands as 
a function of a. In this case a higher coupling strength (Uq = 1.2meV) 
yields important Landau mixing. For a strong magnetic field the electric 
field contribution is small and the spectrum is very similar to that previously 
obtained by Petschel and Geisel [SHI- Three regimes can be identified in the 
plot: (z) In the strong magnetic field limit (small a), the Landau bands are 
well separated and re-scaled butterfly structures are clearly observed, (ii) In 
the intermediate region (er ~ 1) the periodic potential induces Landau level 
overlapping, (m) as a increase very narrow levels develop as a result of the 
electric field, this effect is clearly depicted in the right inset of Fig. |21 

The effects produced by the electric field can be clearly identified in Fig. 
El and Fig. where the energy spectrum as a function of the electric field 
intensity (p = eaE/Uo) are displayed. In the first case (Fig. EJ) Landau 



18 



mixing is neglected and the results for a single band and various selections 
of a are displayed. We first notice that for a = 1/2 and E = the two 
band are degenerated at £ = 0, this degeneracy is lifted for a = 1.0038/2; a 
further increment of a to a — 1.1/2 yields a transformation of the continuous 
bands into a series of quasi discreet levels. For a = 3/2 the band structure 
is observed at small values of E; however as the electric field intensity is 
increased a transition from a continuos band to a discreet level spectrum is 
observed. The corresponding plot in Fig. 0] includes Landau mixing and 
shows the energy spectrum as a function of (p = eaE/Uo) for the first two 
Landau bands. For a — 1/2 every Landau level splits in two bands, the 
degeneracy is lifted by the finite Landau mixing. Again as E increases we 
observe that these bands evolve into a series of quasi discreet levels. In 
these plots we observe that the evolution from the band to discreet structure 
of the spectrum produces a interesting pattern in the distribution of gap 
regions, this structure can be related to a series of avoided and exact crossing 
effects of the quasienergies; these effects become more evident in Fig. Eh- It 
was pointed out in reference jHIj that the coexistence of avoided and exact 
crossing of quasienergies in the Floquet spectrum can be understood on the 
basis of a generalized parity symmetry; in such a way that states belonging 
to the same parity number develop avoided crossing, but those of different 
parity symmetry can yield exact crossing. 

In Fig. El the density of states as a function of the ki pseudomomentum 
and the quasienergy £ is displayed. The density of states was calculated 
from Green's function discontinuity obtained from the continued fraction 
expansion. Three values of o = 1/2,3/2,5/2 are selected, consequently two 
bands appear when E — 0. In the superior panels it is observed that both 
as E increases or B decreases, the spectrum for a fixed value of &2 evolves 
into a set of narrow bands. When &2 is varied the spectrum becomes almost 
continuous, except for the small gaps that open between the bands. 

Based on the previous results it follows that the EMB model preserves 
the band structure. For weak electric fields the bands are grouped forming 
the "butterfly spectrum"; whereas as the intensity of E increases the bands 
are replaced by a series of quasi discreet levels: a "magnetic a Stark ladder" . 

4.3 Adiabatic approximation 

In this section we exhibit an alternative expression for the effective Schrodinger 
equation that governs the dynamics of the system. Based on this formalism 
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we shall find an approximated adiabatic solution, that throws further insight 
into the physical results obtained in the previous section. Let us consider 
the generalized Harper equation ([48)1. taking into account that the matrix 
Q m is periodic in m with period p, we find convenient to define the unitary 
transformation 

= U ^(^) c - > U ^ (0) = J^eWWA 5 ltmmod[p] , (51) 

m 

where [m/p] denotes the integer part of the number, and the delta enforces 
the condition I = m, mod[p\. Notice that whereas the index m in c m runs over 
all the integers, the corresponding label of the state di(<f>) takes the values 
I = 0,1,2, ... ,p — 1. The new vector state di{4>) satisfies the periodicity 
condition 

d l (<f>+—) = (52) 



qa J 

and the transformation matrices fulfill the following properties 

oo 

U I>m fo>) Ulj, (0) = S hl> 8{<t>-<p), 

m=—OD 

/ ^ W U li (^) (^) = <W • (53) 
Applying this transformation to the generalized Harper equation yields 



E M {fa + 0, k 2 ) d {fa + 0, k 2 ) = I £ + Ek 2 - iE^- )d{fa + 0, k 2 ) , (54) 



where the Hamiltonian is reduced to a p x p block form 



E M {fa + 0, k 2 ) 



I Q Q + 

Q- Qi Q + 
Q- Q 2 g+ 



V Q~ 







Q~ \ 




Qp-i J 



(55) 
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here Q = Q (ki + 0, k 2 ), Q m = Qm {h + 4>, k 2 ) and Q + = Q + (k t + 0, k 2 ). 
Taking into account that each block in (J55|) is a L x L matrix determined by 
the Landau indexes, we have obtained that the transformation (j51|) reduces 
the infinite dimensional representation of the Schrodinger equation to a pL x 
pL representation. The price is that a non-local derivative term respect 
to parameter has been added. In the absence of electric field, Eq. ()54|) 
represents a finite dimensional eigenvalue problem, that can be solved by 
direct diagonalization. We observe from Eqs. and that appears to 
be directly related the longitudinal quasimomentum hi, in fact if we redefine 
k\ + — > k\ in Eq. (J53j) yields 

#m (h, k 2 ) d (h, k 2 ) = (e + Ek 2 - iE-^-J d {kt, k 2 ) . (56) 

This expression confirms the form of the Schrodinger equation fK->()D previously 
discussed, in which the "dynamics" is determined by a differential equation 
with respect to the longitudinal quasimomentum. However we now have an 
explicit finite-dimensional matrix representation for the Hamiltonian. 

Let the "instantaneous" eigenstates of Hm be with energies A^ a \ i.e. 



H M h {a) (h, k 2 ) = & {a) (k h k 2 ) (k h k 2 



(57) 



where a labels the band state of the non-perturbed problem. A solution of 
Eq. (JoTlj) can readily be obtained in the adiabatic approximation as follows 



d {a) (ki, k 2 ) 



exp 



, 1 

-l — Ki + I — 

E E 



h {a) (k h k 2 ) , 
(58) 



where the Berry phase 7 (0) is determined from the substitution of the pre- 
vious expression in Eq. 







7 («) (k 1; k 2 ) = i d<\> (0, k 2 ) —hM (0, k 2 ) . 



(59) 



The energy eigenvalue is then determined by the periodicity condition ()52jl . 
hence the change of the phase of the wave function (|58jl must be an integral 
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multiple of 2tt, consequently the spectrum in the adiabatic approximation is 
given as 



S^ (k 2 ) = nEaa+^- J**'** dk x A a (k u k 2 ) + ^^ a) (2rr/qa,k 2 ), n = 0, ±1, 

(60) 

This result deserves some comments. The energy spectrum in the pres- 
ence of electric and magnetic field contains a series of discreet levels separated 
by 

hE 

AS = Eaa = — , (61) 
aB 

where we have used Eq. and restored units. These levels are similar to 
the Wannier levels that appear when an electric field is applied to an electron 
in a periodic potential, the energy separation being proportional to aE. In 
the present case, the band structure parallel to the electric field is replaced 
by a set of discreet steps, this "magnetic Stark ladder" are characterized by 
a separation proportional to the electric field intensity, but inversely propor- 
tional to both the lattice separation and the magnetic field. The existence of 
these levels can be explained by the following argument. In the presence of 
simultaneous E and B fields, the electron travels between lattice points in a 
time r = aB/E. As long as the electron does not tunnel into another band, 
the motion appears as periodic, with frequency u = 2tc/t, corresponding to 
a series of energy levels whose separation AS = ftw coincides with the result 
in Eq. ([61)1 . This magnetic-Stark ladder is combined with the Hofstadter 
spectrum represented by the second term in Eq. ()60|) and the contribution of 
the Berry phase, the competition between these three factors was discussed 
in the numerical analysis of the previous section. We notice that the Berry 
phase is written as the integral of the longitudinal component of the Berry 
connection A a (k) discussed by Kohmoto 



5 Summary 

We have considered the quantum mechanics of a periodic electron system 
driven by electric and magnetic fields. We presented a thoroughly discussion 
of the symmetries of the electric-magnetic Bloch problem. These symme- 
tries were utilized in order to construct the wave function of the system. 
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The space-time representation of the wave function (|33I34|) fulfills the Bloch- 
Floquet theorem with modulation functions that satisfy the generalized Bloch 
conditions given in ([35)1: the phases that appears in these conditions are de- 
termined by the dual-gauge potential related to the symmetry operators (J7J). 

The dynamics of the system is governed by a finite difference equation 
(HSJ), that represents a generalization of the Harper equation, the equation 
includes the following effects: (1) an arbitrary periodic potential, (2) an 
electric field orientated in a commensurable direction of the lattice, and (3) 
the coupling between different Landau levels. The model previously discussed 
by Petschel and Geisel [30J, as well as the Harper equation (JTJ) are recovered 
if the appropriated limit is enforced. 

A detailed numerical analysis based on the solution via a matrix continued 
fractions was carried out, and the properties of the generalized Hofstadter 
spectrum were discussed. In the strong magnetic field limit the Landau bands 
are well separated and re-scaled butterfly structures emerge; in the intermedi- 
ate region the periodic potential induces important Landau level overlapping; 
finally as the electric field intensity is increased a transition from a continuos 
band to a discreet level spectrum is observed. In this "magnetic Stark ladder" 
the equally spaced discreet levels are separated in integer multiple of hE/aB. 
We find that in order to preserve the self-similarity structure of the Butter- 
fly spectrum, besides the strong magnetic field condition B > iih/(ea 2 ), the 
electric field intensity should be restricted according to Eq. (joTIj) . 

An effective equation of motion in k— space was derived (|56|) . In this 
equation the "dynamics" is governed by a differential equation with respect 
to the longitudinal quasimomentum, the explicit matrix representation for 
the Hamiltonian is given in ()55|) . Based on this formalism an approximated 
adiabatic solution for the energy spectrum (|6U|) was obtained. This expres- 
sion includes: (1) the discreet "magnetic Stark ladder" term (JoTj) . (2) and 
averaged Hosftadter contribution, and (3) the Berry phase contribution (foTJj) . 
This Berry phase is written as the integral of longitudinal component of the 
Berry connection discussed by Kohmoto [T2j in relation with the quantiza- 
tion of the Hall conductance. We finally remark that the present formalism 
should set the basis for the study of the Hall conductivity beyond the linear 
response approximation. 

We acknowledge the partial financial support endowed by CONACyT 
through grant No. 42026-F. 
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Figure 1: The energy spectrum inside the lowest Landau level as a function of 
the inverse magnetic flux a. The energy axis is rescaled to e = £ / [f/o e_7rcr ^ 2 ] • 
All values of k\ G [0,27r/g6] are included. The parameters selected are: 
m* = 0.067m e , a = 100 nm, U = 0.5meV, E = 0.5V/cm, and k 2 = vr/2. 
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Figure 2: The energy spectrum for the four lowest Landau bands is plotted as 
a function of magnetic flux a. The coupling strength K = 6 yields important 
Landau mixing. All values of k\ e [0,27r/g6] are included. The parameters 
selected are: m* = 0.067m e , a = 100 nm, Uq = 1.2meV, E = 0.5V/cm, and 
k 2 = 7r/2. The right inset shows the details of the spectrum in a reduced 
region that was originated from the lowest Landau band. Seven Landau 
levels are included in the calculation in order to attain convergence. 
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Figure 3: The energy spectrum for a single Landau level as a function of 
the electric field intensity (p = eaE/Uo) for: (a) a — \, (b) a — L0 ° 38 , (c) 
a = hp, (d) a — |. The energy axis is rescaled to e = £/ [?7oe _7rcr ^] . All 
values of k\ G [0, 2n/qb] are included, and the other parameters are selected 
as: m* = 0.067m e , a = 100 nm, Uq = 0.5meV, and k 2 = ir/2. 
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Figure 4: Energy density plot for the two lowest Landau levels as a function 
of the electric field intensity (p = eaE/Uo) for a = 1/2. All values of k\ G 
[0, 2ir/qb] are included, the other parameters take the same values as in Fig. 

El 
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Figure 5: Energy density plot for the lowest Landau levels as a function of 
the transversal pseudomomentum fc 2 . Three values of a = 1/2,3/2,5/2 are 
selected, so two bands appear when E = 0. Three values of E are considered: 
p = eaE/Uo = 0, 0.002, 0.01. 
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